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PACS 31.10.+Z - Theory of electronic structure, electronic transitions, and chemical binding 
PACS 31 . 15.ae - Electronic structure and bonding characteristics 

Abstract. - We examine the energy spectrum of a charged particle in the presence of a non- 
rotating finite electric dipole. For any value of the dipole moment p above a certain critical value 
Pc an infinite series of bound states arises of which the energy eigenvalues obey an Efimov-like 
geometric scaling law with an accumulation point at zero energy. These properties are largely 
destroyed in a realistic situation when rotations are included. Nevertheless, our analysis of the 
idealised case is of interest because it may possibly be realised using quantum dots as artificial 
atoms. 



, Introduction. — The problem of a charged particle in 
the field of a physical electric dipole serves as the starting 
point for the description of a variety of important physical 
processes, ranging from e.g., the passage of muons through 
4 substance [1], the determination of carrier mobilities and 
^charge trapping states in various condensed matter sys- 
^>tiems [2,3], to the chemical bonding of two atoms [4,5]. 
^ Despite the vast amount of literature on the quantum me- 
5^qhanics of this system, there have been remarkably few 
detailed, numerically exact studies [6,7]. Past efforts have 
focused on the narrow region near criticality where scatter- 
ing states from the continuum are brought down into the 
discrete spectrum [8-11]. The primary impetus for these 
studies is rooted in chemistry where the determination of 
the critical dipole moment (s), along with the binding en- 
ergy of the charged particle, are very important. It is well 
established that a non-rotating, rigid dipole can only bind 
a charged particle if the dipole moment exceeds the critical 
Pc « 1.6249 Dcbyc. 

While part of the reason for the heretofore limited nu- 
merical results can be attributed to the computationally 
intensive nature of the calculations, the more significant 
issue has likely been the absence of additional physical 
motivation to warrant further studies. It is now known 
that the inclusion of rotation drastically alters the ide- 



alised spectrum [12, 13] of the stationary dipole. With the 
advent of quantum dots as artificial atoms, however, there 
is a possibility that the rotational degrees of freedom are 
irrelevant. In that case the idealised model of this paper 
may be of interest to the ultra-cold atoms community, who 
look for such a spectrum at Efimov resonance [14, 15]. 

In this Letter, we present a detailed, numerically ex- 
act investigation for the bound-state energy spectrum of 
a charged particle in the field of a stationary finite dipole. 
This model occurred in the physics literature as early as 
1947 [1]. It was firmly established and later "rediscov- 
ered" [16] that a finite stationary dipole can capture an 
electron if the dipole moment, p, is larger than a critical 
value Pc. Furthermore, for supercritical values of p not 
only one but infinitely many bound state appears in the 
spectrum, with an accumulation point at i5 = 0. 

We focus our attention on the same system: a rigid, sta- 
tionary dipole, hence we exclude the effect of rotation of 
the dipole from our calculation. Nevertheless we have to 
note that rotation may play an important role in influenc- 
ing the spectrum of the electron, especially if the moment 
of inertia, /, of the dipole is small. Garrett [12,13] showed 
that the critical dipole moment, Pc, is larger compared to 
that of a stationary dipole and depends on /. Addition- 
ally, for p > Pc the number of bound states is reduced to 
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a finite value and the infinitude of the number of bound 
states is only recovered if / — > oo. 

Until now the striking similarity between the detailed 
spectrum of our system and the three-body Efimov lev- 
els [17,18] has not been emphasised. Despite the different 
origins of these two effects [19], our numerical and analyt- 
ical results of the bound electronic spectrum for p > pc 
yield two features that resemble an Efimov-spcctrum: 

• a geometric scaling law of consecutive binding ener- 
gies, En+i/En = constant, and 

• an accumulation of the bound states as i? — > 0^. 
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Fig. 1: Schematic representation of the model. The line joining 
the center of the dipole to qe is of length r, and makes an angle 
9 with the z-axis. 



The model. — Our physical model consists of three 
charged particles as schematically illustrated in figure [TJ 
We keep the position of two particles fixed (i.e., Qi and 
Q2) and look for the bound state(s), if such exists, of the 
third particle. Specifically, for the stationary rigid dipole, 
we choose Qi = —Q2 = Q, and the third particle is taken 
to be an electron. Therefore the Hamiltonian, H, for our 
system reads 
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It is convenient to introduce the dimcnsionless dipole mo- 
ment, A, with the definition 
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where p denotes the electric dipole moment, p ~ Qd, and 
d is the size of the physical dipole. Finally, scaling all 
lengths as r^/c?, reduces eq. ([T]) to 



H = -V^-|-A 
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In our simplified picture, A characterises the strength of 
the interaction between the originally unbound electron 
and the physical dipole. If A does not exceed the threshold 
value Ac ~ 1.27863, the third particle cannot be bound. 



The minimum A then represents the critical strength at 
which the first zero-energy bound state appears. Note, 
from eq. that for r ^ d, the potential goes like 

cos(6')/r^ in spherical polar co-ordinates, where r and 9 
are defined in figure [T] For r <^ d, however, the potential 
has no singularity, and goes to zero as r — > 0. At threshold, 
the spatial extent of the wavefunction for the third parti- 
cle becomes much greater than the size of the dipole [20] , 
and the short-distance length scale, d, set by the size of 
the dipole may be ignored. As a result, the same critical 
value Ac is obtained for both the finite and for the point 
dipole [21-24]. For p > pc, however, the point dipole po- 
tential cos{9) jr^ has no lower bound in energy because of 
the singularity at the origin [23], and is unphysical. By 
contrast, our finite-sized dipole potential is well-behaved 
for all r, thereby yielding a meaningful bound state spec- 
trum for p > Pc- Moreover, it is the inverse-square nature 
of this potential for large distances that gives rise to the 
Efimov- like features of the spectrum [25] . 

Numerical implementation — The bound state energies, 
E < 0, and wavefunction, 5*, for the electron are obtained 
from the reduced Schrodinger equation. 
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where = —2mf.Ed^ /It? > is the dimcnsionless energy. 
Equation (4) is separable in prolate spheroidal coordi- 
nates [26], resulting in the following equations [ip £ [0, 27r[, 
C e [l,oo[ and 77 € [-1, 1]) 
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where the wavefunction has been factorised as 'J = 
S{£,)T{r])F{(p). The reduction of the three-dimensional 
Schrodinger equation to a set of three one-dimensional or- 
dinary differential equations (ODEs) requires the appear- 
ance of two separation constants, A and m. Separability, 
as the manifestation of a symmetry, also means here that 
these constants are the eigenvalues of conserved quanti- 
ties [27,28]. If there is no electric dipole, i.e. A = 0, the 
electron is free and the solution of the original physical 
problem is known. In this case the two conserved quanti- 
ties are and L2 with eigenvalues £{£+ 1) and to, where 
^=0,1,. . . and TO =0, ±1, . . . ±£. Since our potential docs 
not possess spherical symmetry for A > 0, the {2i + 1) 
degeneracy of to is absent, although the azimuthal sym- 
metry of the finite dipole still guarantees the conservation 
of L^. Notice that the same is not true for L^, therefore 
£ is no longer a good quantum number. In what follows, 
we adopt the notation established previously in the liter- 
ature [6, 7, 29] and label the energy eigenvalues of H by 
a set of three numbers (n^, rt^, to) where counts the 
zeros of S{£,), while rijj does so for T{'i]). 
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Although the general solution of the above equations 
can be analytically given in terms of double-confluent 
Heun functions [30], this approach docs not offer a fea- 
sible way for deriving the energy eigenvalues explicitly. 
Consequently, we fall back to numerical methods [31,32] 
for finding normalisable solutions to these equations for a 
given value of A and . 

In the first instance, we assume that the separation con- 
stant A can be different in the two equations. In this way 
we obtain two sets of curves, viz. vs. Aiad ^^nd ^angj 
for the radial (|5b[) and for the angular equation ([Sc]). re- 
spectively. Only the latter curve, Aang, depends on the 
dimensionless dipole moment A, since A is present in (|5cp . 
but absent in (j5b|) . The crossing(s) of these two sets of 
curves represents the desired solution to the system of 
ODEs since the separation constant A = ^rad = ^ang 
and become common for the two equations. 

Figure [2] shows the two sets of curves, A^^d and Aang 
for A = 0, which case represents a free electron, there- 
fore is conserved and £ is a good quantum number to 
label the curves with. The fan structure of A^ng is ap- 
parent in the figure and remains intact for non-vanishing 
dipole moments. For < A < Ac, there is no common 
point for the two sets of curves in the > domain, 
thus no bound state exists. On one hand, it can also be 
proven, that lim^2^g+ (Aiad) = ^V^i i-o. all Aiad curves 
approach the value (—1/4). On the other hand as we in- 
crease the dipole moment the Aang curves shift towards 
left (see fig. ^ . Therefore there must be a critical value 
of the dimensionless dipole moment A at which the first 
^ang curve, labelled by (n^ = 0,m = 0), also crosses the 
abscissa at ( — 1/4). However, above this critical value of 
A, ^ang(0,0) crosses infinitely many curves of Aiad! this is 
the mathematical origin of the infinite "tower" of bound 
states if A > Ac. Figure [3] depicts ^ang('^j;, "i) only for 
(n,, = 0,771 = 0). Increasing A further, more and more 
branches of the "fan" are shifted to the left and cross the 
value (—1/4), hence one can deduce a critical value of A 
for each "tower" labelled by (n^, m). The value of the five 
lowest Ac(n^,m) are reported in table [1] As A exceeds a 
given Ac(n,;,m) a new "tower" of infinitely many bound 
states opens, e.g. for A = 16 there will be three sets, 
{nri,m) = (0,0), (0,1) and (1,0) each containing infinitely 
many bound states labelled by =0, 1,2,... 

The values of ^ang and Aj-ad are calculated as eigen- 
values of matrices or from three-term recursions [31,32]. 
Numerically, it is demanding to obtain precise values for 
the crossings, particularly for the rapidly decaying A^ad 
curves near the critical region A ^ A^ . 

Results. — For a given set of (n^, to), if A > Ac(n^, m) 
there are an infinite tower of bound states, irrespective of 
the value of A. This result is in marked contrast to what 
is found in the Efimov effect, where the infinite tower of 
bound states for the three-body system occurs only at 
criticality. We note that the lowest critical dipole mo- 
ment Ac ~ 1.2786 corresponds to an electron with quan- 
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Fig. 2: The angular (thin lines) and radial (thick lines) sep- 
aration constants ([Sbjc) are plotted versus the dimensionless 
energy eigenvalue for A = 0. Each curve of Aang, except the 
ones for m = 0, is doubly degenerate, therefore each branch 
originating from l{i + 1) has {2£ + 1) curves altogether. The 
curves of the radial separation constant are labelled by (n^, 
m — 0). If A < Ac there is no crossing at positive k, and there 
is no bound state. 
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Fig. 3; The angular (thin line) and radial (thick line) sepa- 
ration constants are plotted versus the dimensionless energy 
eigenvalue . The curves of the radial separation constant are 
labelled by (n^, m = 0). Their rapid decrease as energy ap- 
proaches zero is apparent. Moreover, figure shows only {n^ = 0, 
m = 0) angular separation constant for four different A values: 
two below the critical value, one at the critical value (~ 1.2786) 
and one much larger than the critical value. The small vertical 
arrow at Arad = ^ang = — 1/4 captures the position where the 
zero-energy bound states first appear. 

tum numbers (0,0,0), that gives a zero-energy bound state. 
For n^i > 0, it is well-known in the literature that higher 
values of Ac arc needed. Our exact numerical results (see 
Table [T|) confirm these super-critical values [24] with high 
accuracy. 

Figure 5] summarises the central results of this paper. 
Focusing first on the curves with symbols, the main figure 
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Table 1: The five lowest values of Xc{n,^,m) are given. If A > 
Ac (n,,,m) a new set of infinite set of bound states is created. 
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shows the ratio of consecutive bound state energies for 
a variety of dimensionless dipole moments. It is evident 
that at criticahty, A = Ac, the ratio diverges, which is 
expected given that this is the dipole moment for which 
the zero energy state appears. In addition, we see that 
above criticality, the ratio 7, decays exponentially. Similar 
curves are also found for other values of Ac (n^,m), and 
when examined in detail, these numerical findings suggest 
a universal form for 7 near criticality. Indeed, following 
the earlier work of Abramov [29], one can deduce a simple, 
universal expression for the ratio of successive eigenvalues 
valid near the critical region {n^ > 0, and m are fixed) 

i.K + i,n„H^^^ / r^(i) \ 

The 7 values obtained from eq. ^ are also shown in the 
main panel of figure [H Note that very near criticality, 
the analytical and numerical ratios are indistinguishable. 
However, as wc leave the critical region, the two curves 
begin to deviate appreciably, with the result that the ana- 
lytical prediction overestimates 7. The linear relationship 
between the and in the inset to fig. [?] clearly estab- 
lishes the Efimov-like scaling for the bound state energies. 

In fig. m we display the dimensionless energy spectrum 
of the electron in the field of the finite electric dipole for 
A = 20, for which 7 k, 5.38. The accumulation of bound 
states as E' — > 0~ is emphasized by magnifying the spec- 
trum around zero energy, and continues indefinitely. 

Outlook. — We propose an experimentally feasible 
scenario to observe this geometrical scaling. Our sugges- 
tion involves positively and negatively charged quantum 
dots [33]. By assembling two such dots (representing our 
fixed charges) onto an insulator matrix one could design a 
near ideal representation of our model system. This reali- 
sation lacks the full rotational symmetry around the axis 
of the dipole, contrary to our model analysed above, how- 
ever it only restricts the possible values of the m quantum 
number to even integers and should not alter the infinite 
number of bound states and their scaling property. The 
loosely bound states could be detected by standard photo- 
detachment spectroscopy [34]. The additional advantage 
of this scenario would be the possibility of tailoring the 
dipole moment to a favourable value. 
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Fig. 4: The ratios of consecutive energy eigenvalues, 7, are 
plotted for m = and m = 1 in a semi-logarithmic graph. 
The solid thin and dash-dotted lines represent the analytical 
approximation (j6]) while the calculated data are represented 
by -I- (m = 0) and by x (m = 1). The dotted vertical lines 
depict the Ac values (see table [TJ. Inset: dimensionless energy 
eigenvalue function of for different values of A in 

a semi-logarithmic graph. Curves are labelled with increasing 
values of A from left-to-right. Lines are a guide to the eye. 




Fig. 5: Dimensionless energy eigenvalues hP' of the (0,0) band 
are shown for A = 20. Accumulation of energy levels as 
( — K^) 0~ is obvious and also demonstrated by magnify- 
ing the spectrum around the zero energy. The scales on each 
axis are different. 

Conclusion. — We have investigated the energy spec- 
trum of a charged particle's bound states in the field of a 
fixed, finite dipole. Using an exact numerical approach, we 
have confirmed the numerical values of previously known 
critical dipole moments required to bind the charged par- 
ticle in different symmetry states. The new result is that 
for each set of fixed quantum numbers {riri , m) there are 
an infinite number of bound states, labelled by nj, for 
which the consecutive energy eigenvalues are ordered in 
geometrical scaling. We refer to this scaling as Efimov- 
like because geometric scaling of the bound state energy 
levels is a signature of the Efimov effect. Wc also find 
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an accumulation of states near zero energy, as in the Efi- 
mov states. However, our system is not equivalent to Efi- 
mov's original scenario. In the classic Efimov effect, three 
identical (neutral) bosons interact pair-wise via a short- 
range potential, with the result that an infinite scries of 
excited three-body energy levels appear only when at least 
two of the two-body subsystems arc at threshold. In our 
model, wc have a long-ranged electrostatic dipolar interac- 
tion, and the three particles are distinguishable. Further- 
more, the geometric scaling and accumulation of energies, 
i? — » 0~, in our system occurs for any A > A(7i,,, m). Nev- 
ertheless, the common origin for the geometric scaling in 
both cases is the presence of an inverse square potential 
in the description of the bound state energy spectrum of 
the system. We have also suggested a possible experimen- 
tal scenario in which the results of this paper could be 
investigated. 
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